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Abstract 



Wheeler-DeWitt equation is applied to k > Friedmann Robertson Walker metric with 
various types of matter. It is shown that if the Universe ends in the matter dominated 
era (e.g., radiation or pressureless gas) with zero cosmological constant, then the resulting 
Wheeler-DeWitt equation describes a bound state problem. As solutions of a non-degenerate 
bound state system, the eigen-wave functions are real ( Hart le- Hawking) and the usual issue 
associated with the ambiguity in the boundary conditions for the wave functions is resolved. 
Furthermore, as a bound state problem, there exists a quantization condition that relates 
the curvature of the three space with the energy density of the Universe. Incorporating a 
cosmological constant in the early Universe (inflation) is given as a natural explanation for the 
large quantum number associated with our Universe, which resulted from the quantization 
condition. 

It is also shown that if there is a cosmological constant A > in our Universe that 
persists for all time, then the resulting Wheeler-DeWitt equation describes a non-bound 
state system, regardless of the magnitude of the cosmological constant. As a consequence, 
the wave functions are in general complex (Vilenkin) and the initial conditions for wave 
functions are a free parameters not determined by the formalism. 



I. INTRODUCTION 

It is well known that the Hamiltonian formulation of classical general relativity is a con- 
straint system [1,2]. Mathematically, the existence of constraints can be understood as a 
consequence of the fact that not all the components of the metric tensor g^ v are not dynam- 
ical coordinates. In the Arnowitt, Deser, Misner (ADM) decomposition of the space-time 
metric [3], the constraints arise because canonical momenta conjugate to lapse function and 
shift vectors vanish identically. Quantization of a system via Hamiltonian goes under the 
name of canonical quantization, in which one has the option of using either Schrodinger 
or Heisenberg representation of operators. Vanishing of the momentum conjugate to the 
lapse function, in Schrodinger representation, gives the celebrated Wheeler-DeWitt equation 
[4,5]. Here, operators act on "superspace" that is an infinite-dimensional space of all possible 
3-geometries, which are not related by diffeomorphism, and all matter field configurations. 
In practice, to make problems analytically tractable, the infinite dimensional superspace is 
reduced to a finite dimensional mini-superspace by working in models with idealized symme- 
tries [6] . In cosmology, the metric of primary interest is that of Friedmann Robertson Walker 
(FRW) Universe. The gravitational sector of the corresponding mini-superspace models have 
only one degree of freedom, which is the scale factor of the Universe. 

The operator ordering problem, which arises in the formalism, has been addressed in the 
literature [7-9]. Halliwell [7] and Misner [8] have demonstrated that the operator ordering 
ambiguity in the Wheeler-DeWitt equation can be completely fixed by demanding invariance 
under field redefinition of both the three- metric and the lapse function for D > 1, where 
D is the number of gravitational coordinates in a mini-superspace model. In semiclassical 
calculations, operator ordering is usually dictated by convenience. This is because a wrong 
factor ordering introduces an error that is in higher power of Planck's constant. 

A scalar field is usually used in the matter sector of the Wheeler-DeWitt equation [10-13]. 
This is because one is particularly interested in applying the Wheeler-DeWitt equation to 
elucidate properties of our early Universe near the Planck epoch, when classical description 
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of matter and gravity must necessarily break down. 

Because the Wheeler-DeWitt equation is a second order differential equation, a priori 
the initial conditions (or boundary conditions) are not determined by the formalism and 
must be specified to completely describe a system. Hartle and Hawking [14,15] have pro- 
posed "no boundary" boundary condition, for which the resulting wave functions are real. 
Hartle-Hawking condition is essentially a prescription for selecting four-manifolds that will 
be included in the Euclidean path integral over history. On the other hand, Vilenkin [16,17] 
has proposed a "tunneling" boundary condition, for which the resulting wave functions are 
complex. Vilenkin's condition is a proposal that the wave function should consist solely of 
outgoing waves from the initial singularity. 

In this paper, we go against the conventional wisdom and first study the mini-superspace 
model corresponding to closed FRW metric with zero cosmological constant and p = p ra d + 
Pdust- Several comments are in order. First, it is not physically inane to use a classical de- 
scription of matter in a quantum mechanical equation. The situation is analogous to solving 
for atomic wave functions while using a classical coulomb interaction. Second, because of 
the usage of classical description of matter, the resulting wave functions are adequate in 
describing only the late phase of a universe. There are several motivations for studying this 
particular model. First, the formalism for quantum cosmology is still in its rudimentary 
stage, e.g., the role of time, interpretation of wave functions, usage of Euclidean continua- 
tion, etc. are not clear [18]. Since the correct interpretation of the theory should give the 
classical limit, we feel that it is necessary to look into the consequences of the formalism 
in this classical regime. Second, limiting ourselves to a closed FRW metric is for a techni- 
cal reason. Of the three possible FRW metrics, i.e., closed, flat, and open, only the closed 
FRW metric has a finite action. Third, even though this is a classical regime, application of 
quantum mechanical formalism can lead to interesting results. Indeed, because the result- 
ing Wheeler-DeWitt equation resembles a bound state problem, the quantization condition 
gives an interesting relation between the energy density of the Universe and the curvature 
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of the three space. Fourth, if our Universe is a closed universe that has gone through an 
inflationary phase, in its early epoch but with zero cosmological constant today, then we 
may be able to say something about the nature of a wave function immediately following the 
inflationary phase by invoking the continuity of the wave function and its derivative. Here, 
we get the conclusion that the appropriate boundary condition for a wave function is that 
of Hartle and Hawking. The conclusion is based on the fact that eigen-wave functions for a 
nondegenerate bound state system are necessarily real. 

The layout of the paper is as follows. In Sec. II, some notations and basic results are 
quoted from the literature. In Sec. Ill, we solve Wheeler-DeWitt equation in a closed form 
for the closed FRW metric with p = p rad to get the complete set of eigen-wave functions. 
The analysis is carried out for the Laplacian and for an arbitrary operator ordering. In Sec. 
IV, the analysis is repeated for the more general case of p = p ra( j + Pdust- In Sec. V, we study 
the consequences of including a cosmological constant to the Wheeler-DeWitt equation for 
the k > FRW metric with p = p ra d + Pdust- In Sec VI, we briefly discuss generalization of 
including more gravitational degree of freedom and its anticipated effects on the nature of 
wave functions and their corresponding boundary conditions. 

II. PRELIMINARIES 

The starting point is action for the system, 

I = I M + I G . (2.1) 

Im is an action for matter, and 

is the familiar Einstein's action for gravity. We will be primarily interested in the quantum 
mechanics of FRW metric. 
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ds 2 = -dt 2 + a 2 [j^r 2 + r 2 dn 2 3 j , (2.3) 

where (a) is the scale factor of the Universe, dVL 2 is the metric on a unit three-sphere, and 
k>0,k = 0,k<0 correspond to closed, flat, and open geometry, respectively. 
The "time-time" component of Einstein's equation gives 



a +k — -^—pa 2 . (2.4) 
o 

For FRW metric, the Ricci scalar is 



*-^'i+(=)"4)- M 



For the closed FRW metric, the integral over the spatial coordinates in the expression for Iq 
(2.2) is finite and can be explicitly done. 



3V» 



rfta (a 2 - fc) , (2.6) 



where V 3 = 2n 2 is the volume of a three dimensional space with constant positive (k — 1) 
curvature. We have also integrated by parts. An astute reader may question the justification 
for integrating by parts. The reasoning is quite involved [2,19]. In turns out that there are 
two methods of deriving Einstein's equations (e.g. (2.4)) from the principle of least action. 
In the first method, when the action (2.2) is varied, one needs to require that both b'g^ = 
and 5g^ v = at appropriate boundaries in order for all the resulting surface terms to vanish. 
This is because, in general, the Ricci scalar involves terms with g^ v . 

In the second method, one starts with an action that differs from (2.2) by a term that 
is proportional to the trace of the second fundamental form at the boundary. Here, to get 
Einstein's equation, one only needs to require that 5g^ v vanish on all appropriate boundaries 
when the new action is varied. Mathematically, the addition of the new term to the action 
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can be viewed as an alternative way of enforcing all the surface terms to vanish in the 
variation. 

In the Hamiltonian formulation of general relativity, it is the latter action that is preferred 
because of the nature of the variational conditions. For the closed FRW metric, the resulting 
action for gravity is precisely (2.6). So finally, the justification for integrating by part is that 
in actuality we are using the latter action. 

From (2.6), the canonical momentum conjugate to scale factor is 



7T„ 



SI G 



3V, 



3tt 



: aa = 



-aa. 



5a 4nGk 3 / 2 2Gk 3 / 2 
In terms of ir a , the "time-time" component of Einstein's equation (2.4) is 



(2.7) 



3tt 



1 



ka 



8ttG 



-pa 



= 0. 



(2.8) 



We will not need the expression for the Wheeler-DeWitt equation in its most general form. 
To quantize the system, we replace 7i a — ► — i-^ in (2.8) to get 



_ _ p d_ p d_ /3tt\ 2 1_ 
a da a da + \2GJ k 3 



8ttG 



-pa 



il>(a) = 0. 



(2.9) 



The parameter p represents the ambiguity in the ordering of a and 4-. p — 1 is the Laplacian 



ordering [9,20,21]. In the subsequent use of WKB approximation, a convenient p = 2 will be 
used. 

To describe the matter sector, we substitute 



8nG 



8nG 



Pr 



+ 



A 

3' 



(2.10) 



a a 

where p r , pa, and A are radiation energy density, matter (dust) energy density, and cosmo- 
logical constant, respectively. 

It is convenient to define several dimensionless parameters. 



(2.11) 
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8ttG 



,-, , 37T n /8nG | ..... . 



In terms of these variables the Wheeler-DeWitt equation is 



-a- p ^-aP^- + a\a - 7 ) 2 - /ia 4 - (3 2 
oa oa 



We further define rescaled variables. 



a = aa 



P = 13 /a 



7 = «7 



The Wheeler-DeWitt equation now becomes 



^(a) = 0. 



2U{a) = (a - 7) 2 - /i5 4 - /3 2 



III. MODEL WITH k > 0, p = p 



'rad 
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In this section, we would like to study the Wheeler-DeWitt equation for a closed FRW 
metric with just radiation as matter. There are several reasons for studying this simplified 
case. First, it is solvable in a closed form. Second, while solving this model, we will encounter 
general properties of quantizing a bound state system, i.e., existence of a quantization con- 
dition and no need for an arbitrary boundary condition for the wave function. 

The dimensionless parameters (2.16-2.19) are now 7 = p, = 0, (3 7^ 0, and the Wheeler- 
DeWitt equation (2.20,2.21) is 



^ + P ^ 2U(~) 
da 2 a da 



i,{d) = 0, (3.1) 



with (Fig. 1) 



2U(a) = d 2 - f3 2 . (3.2) 

The dominant behavior of ip for a 00 is e^" 2 / 2 . Factoring out the exponential behavior, 
■0(a) = g(a)e~ a2 / 2 , g(d) satisfies 

% + £-2~a& + 2 - 1-P)g = 0. (3.3) 
da A a da 

Power series method reveals that for p ^ the function g(a) is an even function. Therefore, 
in terms of x = a 2 , g(x) satisfies 

d 2 q ,1 + p x dq .8 2 — l — p. . „, 

We first seek solution for the Laplacian factor ordering (e.g. p = 1) [9, 20, 21]. For the 
wave function tp to remain finite for x — > 00, the parameter (5 must satisfy a quantization 
condition 

(3 2 = An + 2 n = 0,1,2... (3.5) 

The "eigenfunctions" are the familiar Laguerre polynomials [22] 
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g(x) = L° P _ 2 {x). (3.6) 

n ~ 4 

And finally, the wave functions are 



Ma) = e- d/2 L°_ p _ 2 Ca 2 ). (3.7) 
In terms of radiation energy density, the quantization condition is 

The eigenfunctions ip n , corresponding to first few n states are plotted in Fig. 2. 

Before leaving this model, we would like to investigate the consequences of quantizing 
the same system for another operator ordering (p = 2). For p = p ra d, the p = 2 ordering is 
also exactly solvable. The comparison of the two sets of solutions will lend credence to the 
conventional wisdom that operator ordering may be dictated by convenience. 

For p — 2, the wave equation (3.1,3.2) can be rewritten in a form that resembles a 
one-dimensional wave equation with a potential. g(a) defined as ip{a) = g(a)/a satisfies 
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g(i) = 0, (3.9) 

which is the familiar Schrodinger equation for a harmonic oscillator. The only difference 
is that among the complete set of Hermite polynomials, which makes up the solution to 
harmonic oscillator problem, only the odd Hermite polynomials are acceptable in this case. 
This is necessary in order for the wave function ip(a) = to remain finite at a — 0. The 
familiar quantization condition for harmonic oscillator gives 



f3 2 = 2(odd) + 1 = 4n + 3, n = 0,1,2... (3.10) 

Comparison with the result for the Laplacian ordering (3.5) shows that the eigenvalue spec- 
trum is essentially unchanged. It can also be shown that for an arbitrary p > 1 operator 
ordering, the quantization condition is 
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(3 2 = 4n + l+p, n = 0,l,2. 
IV. MODEL WITH k > 0, p = p md + p dust 



(3.11) 



In this section, we would like to study the Wheeler-DeWitt equation for a more realistic 
case of FRW metric with a combination of radiation and pressureless dust. 

The dimensionless parameters (2.16-2.19) are now 7 7^ 0, (3 7^ 0, and fx — 0. To keep 
the model as realistic as possible, we will assume that p ra d/pdust ~ O(10~ 4 ) at a Q = 1, which 
implies that /3 ~ 7. 

The Wheeler-DeWitt equation (2.20,2.21) is now 

V>(a) = 0, (4.1) 
with 

2f/(S) = (S- 7 ) 2 -/3 2 . (4.2) 

A schematic graph of the "potential" U(a) is given in Fig. 3. The solution to this 
differential equation doesn't seem to be analytically tractable for the Laplacian operator 
ordering (p = 1). We will therefore resort to solving the equation for another ordering (i.e., 
p = 2) where the problem is analytically tractable. 

Again for p — 2, the wave equation (4.1,4.2) can be rewritten in a form that resembles a 
one-dimensional wave equation with a potential. g{a) defined as ip{a) = g{a)/a satisfies 

9(a) = 0, (4.3) 

which resembles a simple harmonic oscillator with a shift Ax = —7 in the potential. The 
key difference is that here the range of coordinate is a G [0, 00] whereas for the harmonic 
oscillator the range is [—00, 00]. 
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p d 
a da 



A solution to (4.3) which vanishes for very large positive values of a is 

g(a) = Dp_i(a-i), (4.4) 

where D u (z) is the parabolic cylinder function [23]. 

The quantization condition is obtained by requiring that 

g(a = 0) = Dp_i(-nr) = 0. (4.5) 

This is again necessary in order for the wave function ip(a) = g{a)/a to remain finite at 
a — 0. Disappointingly, it is difficult to calculate the parameters /3, and 7 for which (4.5) is 
satisfied. We are left to resort to WKB approximation to find the approximate "eigenvalues" . 
The appropriate WKB phase integral relevant for (4.3) is 

J= [ P+ \-2U(a)) 1/2 da = ir(n + 3/4:) ra = 0,l,2.. (4.6) 
Jo 

The usual 1/2 is replaced by 3/4 because a = is not a usual classical turning point but an 
end of the coordinate, and also because we must demand that g(a — 0) = [24]. For the 
limiting case of ft — 7, which is the statement that p ra d = and p = pdust, the phase integral 
can be readily evaluated to give J = |7 2 . 

In terms of p = Pdust, the quantization condition is 



)(^P^) 2 = 2(n + 3/4). (4.7) 



y 2Gk ,y 6k 

For our Universe, (3 7^ 7 but 



P T _ Pr 2A; 



O(10" 4 ). (4.8) 



7 2 p d H 2 Q +k 

In the Appendix , the phase integral for (3 7^ 7 case is evaluated to give J ~ |/5 2 . In terms 
of p ra d and pdust, the quantization condition is 



10 



W(^ a ° + <i^>)= 2 <" +3 / 4 >' (4 ' 9) 

If we make a rough estimate that k ~ H%, then the quantum number of our Universe is 

N „ {Gk) -i „ ( three ra dius V ^ / V „ 1q122> ( } 

v ; \ Planck length/ \ Planck length/ v 1 

A comment is in order. It is obvious that the entire analysis was facilitated by the use 
of p = 2 operator ordering. It is also easy to convince oneself that even if we could solve 
for the p — 1 operator ordering the two quantization results should not differ by much. The 
argument goes as follows. Because the operator ordering parameter (p) appears with a -1 in 
the Wheeler-DeWitt equation (4.1), a wrong value for p should affect the small range of a the 
most. But for a universe with a combination of radiation and dust, small value for the scale 
factor corresponds to the radiation dominated era. And since the two exact quantization 
conditions for p = 1,2 orderings with p = p rad are essentially the same (3.5,3.10), we can 
also safely conclude that the quantization conditions for p = 1,2 orderings should not differ 
by much even for the present case. 

Interestingly, the unusual largeness of the quantum number of our Universe can be easily 
explained within the framework of inflationary universe models. This will be addressed in 
the next section. 

V. MODEL WITH k > 0, p = p I&d + p dwb , and A ^ 

In this section, we would like to study the effects of including a cosmological constant 
to a FRW universe which in its late stage of evolution becomes radiation and then matter 
dominated. A straightforward method, at least formally, would be to include a spatially 
homogeneous scalar field with some self coupling V((f>), whose non-zero expectation value 
would play the role of a cosmological constant [25]. Moreover, as a reheating mechanism, 
there should also be some coupling between the scalar field and the radiation field, so that 
after the e 60 expansion the scalar field could transfer its vacuum energy to the radiation field. 
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Needless to say this is a very tall order. Here we make a loose-shoe approximation by putting 
in by hand a cosmological constant Aqut- And to simulate the reheating mechanism, we set 
Agut = &wGp rs ,d at some high redshift (e.g., z = zqvt) and then set Aqut = for z < zqvt- 
The Wheeler-DeWitt equation is now (2.16-2.21) 



A schematic plot of U(a) is given is Fig. 4. From the figure, we can readily conclude that 
it is again a bound state problem. Moreover, we can immediately infer an essential difference 
between including and not including the reheating phase. If there were no reheating but a 
constant A ^ for all time, then the potential U (a) would describe a non-bound state prob- 
lem Fig. 5. Moreover, A — > has a unique role in the parameter space. Intuitively, if Agut 
settles to some small but a finite value A fina i, then the metric can re-enter into a de-Sitter 
phase, albeit the two de-Sitter phases may be separated by a classically forbidden region. 
Mathematically, this is because A appears in U (a) with the a 4 factor, which eventually comes 
to dominate U(a). In the Wheeler-DeWitt formalism, this would correspond to a tunneling 
model, which is analogous to a quantum mechanical description of Alpha decay of a heavy 
nucleus. A schematic plot of classically allowed regions in U (a) for this case is given in Fig. 



If we assume that we live in this classical allowed region with a small Afi na i ~ 0(1)H%, 
then it is very interesting to estimate the tunneling probability. The Euclidean action for 
tunneling region is [26] 



-j^ + ?-£-2tf(5) #i) = 0, 
da 1 a da 



(5.1) 



with 




(5.2) 



6. 



S ft ~ — 



(GA final ) 



-i 



-G- l H- 2 w -122. 



(5.3) 
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And the tunneling probability is [27-30] 



P T oc exp(-\S E \) ~ exp(-122). (5.4) 

Now let us return to the case where A G ut 7^ only at large redshifts (e.g. z > ^gut) and 
Afmai — 0. We are interested in its effect on the quantization condition (4.9). It is easy to 
trace the large n ~ 10 122 back to the factor Gk. In light of this fact, having a cosmological 
constant in our early universe can be interpreted as the mechanism which "generated" this 
large quantum number. This is because one of the triumph of an inflationary universe model 
is its ability to naturally explain the flatness of the observed Universe. 

Finally, let us investigate the nature of ip at early Universe. For z <C ^gut, the wave 
function should be that of p = p ra d + pdust model (4.4). For z <C zqvt Wheeler-DeWitt 
equation for the Laplacian ordering is (5.1,5.2) 



-^ + l^ + a 2 - pd A = 0. (5.5) 

da 2 a da 

Moreover, near a ~ the cosmological constant term is negligible compare to the curvature 
term. For this region, the differential equation has a single independent solution 



V>(5) = Io(~a 2 /2). 

I v {x) is a modified Bessel functions [22]. Its power series expansion is 



(5.6) 



2s+v 



(5.7) 



J " ( * )= §ssW(!)' 

Interestingly, ip(a — 0) ^ because Iq(0) = 1. 

There are also inflationary universe models where p ra d ^ before the inflationary phase 
[31,32]. In such a case, the Wheeler-DeWitt equation for z > z GUT is (2.16-2.21) 



d 2 1 d 

-p^ + --7Z + a 2 -pa A -p 2 
da A a da 
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tp(a) = 0. 



(5.8) 



Again the cosmological constant term is negligible. The resulting differential 

equation has a single independent solution (3.7) 

M*) = e- d2/2 L°_ E ^(a 2 )- (5-9) 
And because L n (0) = 1 the wave function ip n {a = 0) ^ 0. 

VI. DISCUSSION 

In this section we would like to briefly discuss an interesting generalization. Consider a 
(k > 0) mini-superspace model, which ends in p = p ra d + pdust, but now the corresponding 
Wheeler-DeWitt Hamiltonian is invariant under some symmetry transformation by having 
more than one gravitational degree of freedom, e.g., homogeneous but anisotropic Kasner 
metric [33-35]. In general, symmetry is a sign of degenerate in eigenvalues. In such a case 
we would have an interesting consequence. First, the problem is still a bound state, so there 
would be no need for an arbitrary boundary condition when solving for the wave functions. 
But from general quantum mechanical considerations, eigen-wave functions for a bound 
state is guaranteed to be real only if there isn't any degeneracy. So we have a combination 
of Hartle-Hawking's no boundary condition with Vilenkin's complex wave functions with it 
non-zero flux, J = ^((j)Vip* — The situation would be analogous to n 2 degeneracy 

for each E n level in a hydrogen atom. 

VII. CONCLUSION 

In this paper, we have studied the solutions to Wheeler-DeWitt equation for k > Fried- 
mann Robertson Walker metric with various types of matter. First, we have shown that if 
the Universe ends in the matter dominated era, then the resulting Wheeler-DeWitt equation 
describes a bound state problem. As solutions of a non-degenerate bound state system, the 
eigen-wave functions are real (Hartle-Hawking) and the usual issue associated with the ambi- 
guity in the boundary conditions for the wave functions is resolved. Furthermore, as a bound 
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state problem, there exists a quantization condition that relates the curvature of the three 
space with the energy density of the Universe. Incorporating a cosmological constant in the 
early Universe (inflation) was given as a natural explanation for the large quantum number 
our Universe, which resulted from the quantization condition. Second, we have shown that 
if there is a cosmological constant A > in our Universe that persists for all time, then 
the resulting Wheeler-DeWitt equation describes a non-bound state system, regardless of 
the magnitude of the cosmological constant. As a consequence, the wave functions are in 
general complex (Vilenkin) and the initial conditions for wave functions are free parameters 
not determined by the formalism. 

APPENDIX 

We would like to show that 



J = J pdx — (3 



with 

P={p 2 -(l-~a) 2 ) 1 ' 2 . 

For our Universe 

^"^—-0(10-% 

7 2 # + k Pdust 
Therefore, defining x — a/7, we have 



J = J 7 (p 2 - (7-a) 2 ) 1/2 dx 

J i \ 7 7 / 7 
~ fj^lf-J^ dx 
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7 2 



f 



1 (^) 2 -l 

(l_ x 2 ) l/2 + A V 



2(1 -rr 2 ) 1 ^ 



7T/2 + 7T/2 f(|) 2 "l N 



-7T/3 2 . 
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FIGURE CAPTIONS 

1. Graph of U(a) for p = p rad and k > 0. 

2. Graphs of ip n (a) for first few quantum states. Figures 2(a)-2(d) correspond to n — 
[0, 1, 2, 20], respectively. The normalization is arbitrary. 

3. Graph of U(a) for p = p ra d + pdust and k > 0. 

4. Graph of U(a) for p = p rad + p dust , fc > 0, and A GUT ^ for z > z GVT . 

5. Graph of U(a) for > and A 7^ for all time. 

6. Graph of U(a) for p = p rad + p dust , k > 0, A G ut, and A fina i. A G ut is responsible for 
inflation and Ag na i is some final value the cosmological constant has settled into in a 
universe. 
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